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Abstract. 

Difficulties in obtaining good phantoms, improvements in technologies of voxel 
localization, better sequences for water and fat suppression has brought us to 
define a minimal Protocol of home-made quality controls of MRS systems. 
Measurements, defined in the proposed protocol, have, as main goal, to establish 
if peaks quantification predicts realistic concentration values, meaning that, the 
occurrence of this event is a sufficient condition to declare that MRS system 
works good. Moreover, stability measurements helps in a correct data 
understanding. It is, indeed, realistic to think that environmental condition can 
introduce casual errors in the working good system. Discrepancies in the working 
good condition, under stochastic variability (environment), have to be related to 
systematic errors introduced by the set of pre and/or post-processing operations 
and/or by any forms of MRS bad-working tool that differs from the previous. 
The quality control minimal protocol has been executed on a Philips-Achieva 
MRS system utilizing a phantom supphed by the manufacturer. The minimal 
protocol consists of two steps: reproducibility and performance tests. The 
reproducibility of the MRS measurements helps in quantifying the stability of the 
system. The performance test enables to establish if the system is able to 
reproduce concentrations in a realistic way. In both cases good results have been 
obtained: fluctuations of measured values are below 9% and quantification of 
concentration is consistent with the known values. 
1. Introduction 

Magnetic resonance spectroscopy is a powerful diagnostic tool for gaining in-vivo biochemical 
information and has been progressively developed over the last twenty years. MRS techniques, 
combined with standard magnetic resonance locahzation method, can now be used on clinical 
whole body scanners to classify and quantify the metabolites characterizing organic tissues [1]. In 
particular, MRS aims to distinguish among normal and pathological conditions and anables to 
improve non-invasively a diagnostic tool in cancer detection [2,3,4,5,6]. Good implementation of 
MRS acquisition procedures and post-processing analysis depends on the equipment used for the 
experiments. Physicists are, thus, called to perform MRS pre-in-vivo test in order to detect, resolve 
and prevent problems of poor-quality performance, thus enabling the highest possible standards to 
be maintained [7,8,9,10]. The quality tests are necessary to keep account for the in-vivo plaimed 



use of the equipment. In fact, depending on the in-vivo exams their specific acquisition sequences 
have to be tested in phantom study. This approach does not prevent the possibihty of errors 
originating and, specifically, characterizing in-vivo environment but offers a good test of the 
selected sequences and of the quantification performance to give the possibility to exclude the 
fonts of errors in the successive in-vivo inspections. In this paper we propose a minimal Protocol 
of quality controls for MRS systems. The measurements defined in the proposed protocol have, as 
main goal, to establish if peaks quantification predicts realistic concentration values, meaning that, 
the occurrence of this event is a sufficient condition to declare that MRS system works well. 
Moreover, stability measurements help in a correct data understanding. It is, indeed, realistic to 
think that environmental conditions can introduce casual errors (se ci sono errori casuali il sistema 
non pud funzionare bene). Discrepancies in the working well conditions, under stochastic 
variability (environment)(se non ci sono errori casuali=stocacisti=imprevedibili, gli errori vanno 
addotti al pre-processing malfatto), have to be related to systematic errors introduced by the set of 
pre and/or post-processing operations and/or by any forms of MRS tool bad-working that differs 
from the previous one. In the first case error can be identified and solved. The minimal 
measurement protocol has to be repeated in order to confirm that corrections in post and/or pre- 
processing operations have had a positive outcome. In the second case error has to be searched in 
voxel localization, in contamination of the signal by the boundary of the voxel and quality controls 
have to be performed in a deeper way [7,8,9,10]. 



2. Theoretical Background 



In proton nuclear magnetic spectroscopy (IH-MRS) the intensity of the signal is proportional to 
the volume of interest and to the number of nuclei generating the signal (questo non e vero in 
generale, vedi spettro al carbonio). This statement, translated in Irequency domain, means that the 
area under a specific peak is proportional to the number of nuclei (prtons, sono I protoni nucleari a 
precedere) precessing at that frequency. The analytical equation describing this dependence can be 
shown as [1]: 



A~NV 

where A indicates the area under the pick of the signal, V the volume of acquisition and N is the 
number of nuclei (protons) of a specific metabolite. This relationship is the theoretical basis for 
metabolites quantification. Unfortunately, the classification and quantification of the signal 
produced by MR system is a complex and rather technical issue: imperfection in acquisition and 
post-processing tools can introduce deviations from the theoretical behavior. 
To test the linear relationship that links the number of protons with the areas of the peaks we 
decided to monitor the relative ratio among the areas of different known metabolites that are inside 
the phantom. Given two metabolites characterized by a number of protons A/, and Nj^ amplitudes 
signal area A, and Aj in the same volume of interest, the dependence between areas and 
concentrations imphes that: 

— = — (1) 

The knowledge of the concentrations of the phantom allows the experimental assessment of the 
previous equation giving us the possibility to verify the quality of our instruments. We, 
focusedour attention on the signals of water and acetate because they are not J coupled and 
therefore simpler to quantify. 

The measured area depends on the choice of TE and TR. Under ideal conditions (TE = and TR = 
oo) the transverse magnetization decreases and the longitudinal magnetization is fully recovered. 
Otherwise, we must take into account the effects of relaxation and saturation. The relationship that 
links the parameters defined for a PRESS sequence is [1,7]: 
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We named 7\ and 7*2 the longitudinal and transverse decay time and Af-gy and Ayy^i^ the coiTect 
(che e quella vera) and measured areas, substituting Ai and Aj eq. (1) becomes: 
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where Tu and Ty are the 7\ decay times of the two metabolites and ^i^d are their T2 
decay times. 

From equation (2) it is easy to deduce that T2 can be determined by carrying out measures of 
immediate ./Imis a function of TE for a fixed value of TR. Similarly, 7\ can be determined by 
carrying out measures of Aj^^i^ as a function of TR for a fixed value of TE (si che lo e: se tieni 
costante il TR la tua variabile x e il TE. Al variare di x, Amis(=la tua y) varia seguendo una legge 
puramente esponenziale. In particolare un esponenziale smorzato (segno - dell'esponente)). 
Dependence on TE is particularly easy to determine as it is the purely exponential decay. 
Dependence on TR is slightly more complicated and has required the empirical determination of 
the factor 



^0 =^c 

by which eq.(2) can be written in the form: 

-^0 ~ ^mis 



that is in the form of a simple exponential decay. 
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3. Materials 

The phantom used for the measurements is a polyethylene sphere filled with a fluid of known 
chemical composition. 

The spherical phantom contains: 1000 ml of demineralized water, 5 ml of acetate (CH3C00H), 
10 ml of ethanol (CH3 CH20H), 8 ml of phosphoric acid (H3P04), 1 ml of Arquad solution (1%), 
120 mg/ml CuS04. The proton spectrum of this phantom should have the form depicted in the 
figure 1 
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Figura 1 - Previsional spectrum for phantom 



The signal from ethanol CH2 shows a scalar coupling with constant J~ 7Hz. The relative peak has 
a positive sign with a TE = 288 ms and it is inverted for TE = 144 ms as showen in the figure. 



4. Methods 
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The proton spectroscopic experiments 
were performed on a Philips Achieva 
^ I (Philips Medical system, Eindoven, The 

Netherlands) 1.5 Tesla clinical magnetic 
resonance system . As in our Institution 

MRS will be, principally, used for the 

study of breast cancer, the quality control 

1 I test, here described, has been executed on 

the same coil used in in-vivo examination: 

a 7 Channels Phased Array breast coil. The 

spectroscopic phantom was positioned 

inside the right hole of the coil and locked 

i': ; : with two inert plastic compression plates 

to ensure immobility. The pre-scan 

operation has included automated 
shimming and (si local shimming) water 
suppression. [3,4,5,6,11,12] we used a 
single voxel PRESS sequence with water suppression, characterized by TE=288 and TR=1500, to 
reproduce the spectrum of the phantom. A RESS sequence has three slice -selective radio- 
frequency pulses with orthogonal magnetic field gradients and the intersection of the slices defines 
the volume of interest (VOI). The Philips post-processing software for spectroscopy enables to 
process the FID as prescribed by AAPM Report [7]: 



Figura 2 - PRESS 



a) residual water suppression, 

b) phase coiTCCtion, 

c) zero-filling, 

d) apodization filtering, 

e) Fourier transformation, 

f) frequency shift correction 

g) baseline correction. 



In order to quantify concentration, the parameter to be investigated is the area under the peak in 
the frequency domain. Under ideal conditions it coincides with the absorption mode. Whenever the 
initial phase 9 of a FID is not zero, the real and imaginary channels after Fourier transform will 
contain mixtures of absorption mode and dispersion mode spectra: 

i?e(/) = cos 9 A(f) + sin 9 D(f) 
lm(f ) = cos 9D{f )+ sin 9 A(f), 

where Re(^f^, Im,(^f^, A(^f^ e D{f^ are the real, imaginary, absorption and dispersion mode of 

the spectrum and / is the frequency. Phasing a spectrum sorts the real and imaginary channels into 
absorption mode and dispersion mode spectra. The Philips post-processing software supports this 
function in automated or manual ways. The manual phasing, being operator dependent, can 
introduce non-predictable errors. Therefore, we have decided to use automated phasing that, in our 
opinion, gives good results (la parte reale coincide con I'assorbimento e la parte immaginaria con 
la dispersion. Teoria delle Trasformate di Fourier.) 

[M/)=^(/)e/m(/) = £)(/))]. 

Moreover, we have seen that the water suppression sequence combined with PRESS does not 
completely eliminate the water signal. Post-processing with a convolution difference filter can be 
used to ehminate any residual water signal. This filter applies a low-pass filter to the FID, then 
subtracts the filtered signal from the original data in the time domain. 



Finally, spectral resolution has been improved via zero filling: the number of digital points has 
been increased froml024 to 2048 by adding data points of zero amplitude at the end of the signal 
(ok. Scrivo tutto quello che manca). 



5. Results 

4.1 Stability Test 

During a month we have acquired and processed in this way 14 signals to establish their 
fluctuation and the correlated stability of the machine. The average peak area and their standard 
deviation for acetate [A^is {ACE)], ethanoKCHj) [A^is (£"^1)], ethanoKCHj) [A^i^ (ET2)] 
and unsuppressed total water' [^mis (^2^)] 

Amis (ACE) = 1,8 ±0,1 ^mtsC^zO) = 1029,9 ±50,6 

Arr^is (ETl) = 1,7 + 0,2 Arr^is {ET2) = 1,4 + 0,1 

Fluctuation rate of acetate signal is 6% , of water signal is 5%, of ethanol(CH2) signal is 9%, of 
ethanol(CH3) signal is 7%. These values can be consider excellent [7]. The quality control tests on 
long range (a tempi lunghi. E una espressione nota in letteratura.) stability will be executed weekly 
and one will be considered acceptable the measure having a tolerance of two standard deviations 
from the mean value of every metabolite (fitting. Grazie ad un programma chiamato j-miur e 
prodotto da un progetto Marie Curie. Aggiungero dettagli e referenze.) . 

4.2 Quantification Performance 

As already highlighted, the area of the real part of the peak may not coincide with the absorption 
mode. In order to solve this problem, when we tested for stability, we adopted the technique of 
automatic phasing correction, but for the calculation of the decay times, in order to quantify 
concentrations, we have to operate in a different way. We note, first, that the dependence of the 
signal irom damping (lo smorzamento delle diverse componenti di Fourier di uno stesso segnale: 
assorbimento dispersione, parte immaginaria e reale si smorzano con gh stessi tempi Tie T2) is 
the same for the different components of the spectrum and does not depend on the phase. This fact 
implies that the module of the signal decays in the same way. Moreover, the module of the signal, 
by definition, does not contain information on the phase factor. Therefore, the integration of the 
module is equivalent to eliminate all errors related to phase factor. We stress that the quantity 
obtained in this way does not give information on metabolites concentration but allows us to have 
a more precise estimation of the decay times. The plot of the obtained experimental points is 

TE 

shown in figure 3 and 4 (ho ottenuto empiricamente il fattore Aq = A^-gj. e ^2 . L'avevo 
dichiarato nel theoretical backgroimd). 



The post processing Philips software supports a function able to resume the unsuppressed total 
water signal from the suppressed one so that other acquisition result unnecessary. 



T2 measurement: TR=1800ms 
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Figure 3 - y = A^is,x=TE, TE values are: 1000, 950, 900, 800, 700, 600, 500, 400, 300, 200, 
100 ms 



Ti measurement: TE=288ms 
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Figura 4 - y 
1800, 1900, 2000 ms 



X = TR, TR values are: 1100, 1200, 1300, 1400, 1500, 1600, 1700, 



The corresponding and Tj values obtained by experimental data fitting are: 
T2(ACE)=1011,4ms T2(H20)=795,lms 
Ti(ACE)=5398,9ms T i(H20)=3429,5ms 



The calculation of Ti and T2 allows us to find ^cor relative to the signals of acetate and water. 
One must remember, however, that the analysis was developed using area data relative to the 
module of the signal. While, in order to find the ratio among the concentrations we need to 
calculate the area of the absorption components. To do this we use the average values obtained in 
the previous section. 
Therefore, remembering that: 
A^isi^CE) = 1,8 ± 0,1 
and that 



AmisiHiO) = 1029,9 ± 50,6 
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we obtain: 



N(ACE) Aror (ACE) 
V J ^ _^ory J ^ Q QQ24 + 0,0002 

This value must be compared with the ratio between the number N of protons contributing to the 
signal of acetate and water, respectively. Named n{ACE) and n(//20) the number of protons 
per molecule responsible for the signal, m{ACE) and m{H2 0) the molecular masses of the two 
compounds, p(ACE) and p{H20) their density, V{ACE) and V{H20) the volumes of the 
compounds in the analyzed solution, we can thus show; 



N(,ACE) ^ n(ACE)miH20)piACE)V(ACE) 
N(H20) n{H20)m{ACE)p{H20)V{H20~) 

0.002625 

where n{ACE) = 3, n{H20) = 2, m{ACE) = 60, m{H20) = 18, p{ACE) = 

l,0Sgr/cm^,p(H2O) = 0,9gr/cm^,V(ACE) = 5m/andK(//20) = 1000m/. 
The accordance among theoretical and measured value is acceptable: the theoretical value falls in 
the limit of the tolerance interval coupled to the measured one. 



6. Discussion and conclusions 



In the phantom study, here, the coefficient of variance of the detected intensity for signals from the 
Philips standard phantom were lower than 9% within a 1 -month period. The Tj value of water and 
acetate was determined to be 3429,5 and 5398,9 ms, respectively, and the T2 value of water and 
acetate was determined to be 795,1 and 1011,4 ms, respectively. The concentration ratio of water 
and acetate in the standard phantom was determined to be 0,0024±0,0002, which is consistent with 
the value 0,0026 derived from the known volumes of metabolites provided by the manufacturer. 
From the previous analysis we can consider acceptable the quality of data obtained with magnetic 
resonance spectroscopy exams performed on the 16 Channel, Philips Achieva, l.STesla, Magnetic 
Resonance system here analyzed. 

The quality confrol minimal protocol has been executed on a Philips-Achieva MRS system 
utilizing a phantom supplied by the manufacturer. As already highlighted, the minimal protocol 
consists of two steps; reproducibility and performance tests. The reproducibility of the MRS 
measurements helps in quantifying the stability of the system. The performance test enables to 
establish if the system is able to reproduce concenfrations in a realistic way. In both cases good 
results have been obtained; fluctuations of measured values are below 9% and quantification of 
concenttation is consistent with the known values. 



The capacity to reproduce realistic values of concenfrations is, as already highlighted, sufficient 
condition to affirm that the MRS tool coupled to the magnetic resonance image system works well. 
The good working condition of the system here analyzed ensures that the post and pre-processing 
operation have been conducted in a correct way. Other stability measurements have to be 
performed on our machine in order to ensure long term reproducibility of spectroscopic exams (sto 
testando la capacita del sistema di riprodurre valori di concenfrazione= sto testando la 
Quantificazione, unico parametro utile in spettroscopia). 
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